function [SE_DT,SE_DNT]=power_slope(n,C,tau,sigma,pi,f,j,k)

%Initialize distribution statistics
mu_k=zeros(1,length(pi));
s_k=zeros(1,length(pi));

%Calculate distribution statistics
for i=1:length(pi);
mu_k(i)=pi(i) *f(i);
s_k(i)=(1-pi(i))*f(i);
end;

%Standard Errors
SE_DT=(1/(pi(k)-pi(j)))*(n*tau*(1/f(j)+1/f(k))+sigma*(1/mu_k(j)+1/mu_k(k))).^0.5/(n*C)^0.5;
SE_DNT=(1/(pi(k)-pi(j)))*(n*tau*(1/f(j)+1/f(k))+sigma*(1/s_k(j)+1/s_k(k))).^0.5/(n*C)^0.5;

end
